Condensation energy and the mechanism of superconductivity 
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Condensation energy in a superconductor cannot be precisely defined if mean-field theory fails 
to hold. This implies that in the case of high temperature superconductors, discussions of quan- 
titative measures of condensation energy must be scrutinized carefully, because the normal state 
is anomalous and the applicability of a mean-field description can be questioned. A related issue 
discussed here is the precise meaning of a superconducting transition driven by kinetic as opposed 
to one driven by potential energy; we argue that this is a semantic question. 



Introduction: In an earlier paper (CKA) [fij, we raised 
the issue that the notion of superconducting condensa- 
tion energy is ill-defined if the transition cannot be 
described by BCS mean-field theory, where the system 
turns into a normal Fermi liquid with no pairing corre- 
lations once the superconducting order parameter van- 
ishes. Another purpose of that paper was to elucidate 
the interlayer tunneling theory (ILT) j3|. In particular, 
we examined the strong version of ILT proposed by An- 
derson [|], |5|, in which the entire "condensation energy" 
was ascribed to ILT. This proposal turned out to be at 
variance with the c-axis penetration depth measurements 
of Moler et al. [§ in T12201 and was thus falsified. Nev- 
ertheless, we were interested in understanding if it is at 
all possible that ILT plays an important role in enhanc- 
ing the transition temperature, T c , by increasing the bare 
superfluid density, thus defining ILT in a weaker sense, as 
an enhancement mechanism over and above an in-plane 
pairing mechanism JtJ . 

CKA also noted that nominally optimally doped 
T12201 has a specific heat peak [|| that could be approx- 
imately fitted by a two-dimensional (2D) Gaussian fluc- 
tuation contribution to the free energy. This observation 
reflects once again the importance of in-plane pairing cor- 
relations and was an important conclusion of CKA. We 
then asked if there was a sensible procedure to subtract 
the 2D fluctuations and use the remainder of the free en- 
ergy to understand the effect of ILT in the weaker sense 
of enhancement of the bare superfluid stiffness. This was 
difficult, as the correctness and the precision of the spe- 
cific heat measurements |^] were unknown and still are 
because the measurements are yet to be reproduced by 
a second group. In addition, it was not clear over what 
range of temperatures the fluctuation contributions must 
be fitted. Of course, the very notion of Gaussian fluc- 
tuations in a 2D superconducting transition cannot be 
meaningful close to the transition. Despite these difficul- 
ties, an approximate subtraction procedure was used by 
CKA. The result was that the enhancement of the bare 
superfluid stiffness in T12201 was indeed extremely small. 
Nonetheless, we believe that it is conceptually important 
to perform such subtractions, preferably more accurate 
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FIG. 1: Transition temperature across a homologous series: 
HgBa2Ca n _iCu n 02n+2+<5; adapted from Ref. Il2| 



ones, to estimate the effects of ILT. This is also true for 
multilayer cuprates, where superficially ILT seems to be 
important 0, at least in the weaker sense defined earlier. 

To this day, the cause of a striking systematic rise 
and a subsequent drop in T c for a homologous series 
as a function of the number of layers in the unit cell 
is not known. Even if we ascribe the rise to the en- 
hancement due to ILT jl(| , the drop must be ascribed to 
a competing mechanism that develops with the increase 
in the number of layers, perhaps because the inner lay- 
ers have a tendency to become underdoped. A homol- 
ogous series of cuprate superconductors is a family in 
which each member has the same charge- reservoir block, 
but n Cu02-planes in the infinite-layer block, which con- 
sists of (n — 1) bare cation planes and n-Cu02-plancs 
PJ. Clear systematics of T c is only evident within a 
given homologous series. A well studied example E^ ] 
is the family HgBa2Ca„_iCu„02ri+2+(5 whose T c , opti- 
mized with respect to oxygen concentration, as a func- 
tion of n, is shown in Fig. |l|. The formal copper valence 
vcu = 2(n + S)/n is also a bell shaped curve that peaks 
at n = 3. Similar results arc known for other families, for 
which the transition temperatures follow a similar pat- 
tern, often peaking at n — 3 or 4. The issues of the 
dependence of T c on the number of layers and the role of 
ILT remain unresolved. 
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We now address the basic question raised by CKA, 
namely, is condensation energy a precise quantitative 
concept for high temperature superconductors? Given 
the interest that this subject still elicits |)[ [D| [l4|, [If], 
PI [l7| we have decided to publish this brief note 
to elaborate further on condensation energy and on a 
related theoretical question: can the mechanism of su- 
perconductivity be usefully said to be driven by kinetic 
as opposed to potential energy? Our conclusions will be 
that while it is important to identify the mechanism by 
which the condensate is formed, it is a semantic issue as 
to whether or not we describe the transition as driven by 
potential or kinetic energy. 

Condensation energy: Colloquially, the condensation 
energy is the difference of the ground state energies be- 
tween the normal state and the superconducting state. A 
little thought reveals several related problems: (1) What 
do we mean by the normal state? In particular, what if 
there are other broken symmetries Jl^, ||(], ^ll ^2| in 
the regime in which there is no superconductivity and a 
further transition to the unbroken symmetry state at a 
temperature above the superconducting T c ? (2) What if 
the normal state contains superconducting fluctuations? 
(3) What if the normal state changes as a function of 
the magnetic field, or other tuning parameters used to 
destroy the superconducting state? (4) What if the tran- 
sition to the normal state is not a first order transition, 
such that one cannot meaningfully define a notion of a 
metastable state that can be accessed in experiments? 
(5) How should one correctly extrapolate the non-zero 
temperature measurements to T = to access the hy- 
pothetical normal state with the same set of parameters 
for which Nature actually provides us with the supercon- 
ducting state? There are indeed simplifying situations, 
where the complexities mentioned above do not arise in 
the practical sense (2). Thus, when mean field theory 
holds and the normal state is a Fermi liquid with no 
measurable trace of pairing correlations, the simplest ex- 
trapolation of the normal state below T c with the specific 
heat C(T) = jT, where 7 is a constant, is plausible, as- 
suming that there are no other instabilities of the Fermi 
liquid at temperatures below T c . One may further con- 
strain this extrapolation by entropy conservation because 
the difference of entropies between the normal state and 
the superconducting state is zero at the mean field T c 
and at T = @. 

For high temperature superconductors, there are great 
many complexities. The presence of a pseudogap, quite 
unlike a BCS superconductor, makes the extrapolation of 
the normal state (in fact, even its definition) exceedingly 
problematic. If the magnetic field, H, is used as a tuning 
parameter to destroy superconductivity, its large magni- 
tude, H > H C 2, may stabilize some other ordered state 
. Moreover, unlike conventional superconductors, for 
which the effect of the magnetic field in the normal metal 
is a weak Landau diamagnetism, the normal state of high 



temperature superconductors may not be so impervious 
to such high fields necessary to destroy superconductiv- 
ity. The attempt to destroy superconductivity by doping 
Zn to replace Cu suffers from similar problems. In fact, it 
is empirically known that Zn impurities introduce mag- 
netic order in high temperature superconductors p5[ . 

There are even more fundamental reasons for doubt- 
ing the notion of condensation energy. If the transition 
is a continuous transition, there is no way that one phase 
can be continued into the other beyond the transition. 
Therefore, the hypothetical normal state cannot exist for 
the same set of parameters for which the superconducting 
state is more stable; the notion of a metastable state is 
thus not meaningful for a continuous transition. An ex- 
actly solved model illustrates this point beautifully. Con- 
sider the 2D Ising model for which Onsager's result for 
the free energy is known for all temperatures. The an- 
alytic continuation of the free energy, / + , from above 
the ferromagnetic transition point Tp to below Tp was 
obtained exactly by Majumdar ^6|. One gets, close to 
T F , 

/+ — - jr^l (u - up) 2 [In \u - U F \ + in] , (1) 

<±TTUp 

where u = exp(— AJ/ksT), up is its value at the transi- 
tion point, Tp, and J > is the ferromagnetic exchange 
constant. It is seen that the analytic continuation ac- 
quires an imaginary part, which has no obvious physical 
meaning. This is true for any continuous transition for 
which specific heat exhibits a nonanalytic critical singu- 
larity, reflecting a branch point in the complex plane. 
It is even true for infinite order transitions, as in a six- 
vertex model. The exact analytic continuation of the free 
energy of the six-vertex model was obtained by Glasser 
et al. |2?[| . If the transition were instead a first order 
transition, the imaginary part of the free energy could 
be interpreted as the decay of the metastable state [^] . 

It might be tempting to define condensation energy as 
the difference between the exact ground state energy with 
zero order parameter (unbroken symmetry state) and the 
exact ground state energy with a prescribed finite value 
of the order parameter (broken symmetry state). For a 
broken symmetry with a non-conserved order parameter, 
as in a superconductor, this is impossible, simply because 
the order parameter and the hamiltonian cannot be si- 
multaneously diagonalized. To understand the nature of 
the broken symmetry state with a non-conserved order 
parameter, consider the simplest such case: an antiferro- 
magnet for which the staggered order parameter is not 
conserved. In a bipartite lattice, where the Marshall sign 
condition [^9| holds, the ground state is always a sin- 
glet. In a finite volume, the symmetry cannot be broken, 
and, for a large system, the order parameter will precess 
slowly so that no orientation is preferred. The effective 
hamiltonian, 7i c s, that describes this precession depends 
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on the total spin, Stot, and is that of a rotor, given by 

n cS = ^-s? ot - ^s(s + 1) (2) 

where \ — is the total spin susceptibility, in units 

of g/J,B, h = lj XT' is the susceptibility per spin with 
respect to a local uniform magnetic field oriented per- 
pendicular to the staggered order parameter. One can 
imagine deriving this hamiltonian by a renormalization 
group analysis, as the relevant states are all below the 
one-magnon state of the smallest non-zero momentum in 
a box. Even though the actual eigenstates are those of 
total spin, as N — * oo, a tower of excited states collapses 
to the singlet ground state corresponding to S = 0, and 
becomes degenerate with it in the thermodynamic limit 
|30|. The broken symmetry state with a fixed direction 
of the staggered order parameter is a coherent superpo- 
sition in this quasi-degenerate manifold. Thus, the ener- 
getic difference with the singlet ground state vanishes in 
the limit N — > oo. The energetic difference between the 
normal state and the condensed state is identically zero. 

Intuitively, one feels that one should be able to define 
condensation energy variationally. Consider two varia- 
tional wave functions, one of which corresponds to the 
superconducting state with broken U(l) gauge symme- 
try, and the other corresponding to the normal state. Of 
course we have to define what we mean by the normal 
state-a Fermi liquid, a state with another broken sym- 
metry, etc. Similarly, we must also define the order pa- 
rameter symmetry in the superconducting state. Given a 
hamiltonian, we can now calculate the expectation value 
with respect to these states and find the difference in en- 
ergy, hence condensation energy. This is not only model 
dependent, but also calculation dependent. More impor- 
tantly, there is no known experimental method to check 
the correctness of this definition of the condensation en- 
ergy. 

There is one instance in which the condensation energy 
can be defined with little ambiguity [Q, and that is for 
a type I superconductor. In this case, the transition to 
the normal state as a function of a magnetic field is a 
first order phase transition with only a finite correlation 
length. If the normal state is relatively insensitive to the 
applied magnetic field necessary to destroy superconduc- 
tivity, the measurement of the thermodynamic critical 
field H c , as T — » 0, immediately gives the condensation 
energy from the formula 

H 2 

SIT 

where G is the Gibbs free energy, and fl is the volume of 
the sample. Unfortunately, this is unusable for high T c 
superconductors because they are of type II. 

Frustrated kinetic energy: An idea that has been dis- 
cussed often is that the superconductivity in the cuprates 



is driven by the saving of the electronic kinetic energy in 
the superconducting stat e flj ||, [?|, |3l], |||, Q . There are 
some experiments |)| [34| |35j that could be interpreted in 
this manner. 

Cuprates are complex materials with intricate elec- 
tronic structure. If we assume that electron-phonon in- 
teractions do not play a major role, the problem is en- 
tirely electronic in nature. For concreteness, let us as- 
sume that a single band two-dimensional Hubbard model 
is a good effective hamiltonian to understand the low en- 
ergy properties, including the superconductivity of these 
materials. Even if the electronic hamiltonian were more 
complicated, it would make no difference to our basic ar- 
gument. For example, we could also incorporate electron- 
phonon interaction at the expense of making the discus- 
sion more complex. The one-band Hubbard model de- 
scribes processes smaller than energy U and is 

H eS = -t ^ ( C L C 3<r + h - C ') + U n »T"a ( 4 ) 
(ij) * 

The higher energy processes are assumed to be adiabat- 
ically decoupled from the lower energy processes. Here 
Cia- is an electron destruction operator of spin a, and 
is the corresponding density operator. 

When U is large, the model can be reduced to the 
effective hamiltonian called the t-J model, which is 

H t -j= - ^(4c JCT +h. c.) 

m 




with J = 4t 2 /U, together with the constraint < 1. 
The operators c%& still satisfy the fermion anticommu- 
tation rule, but one must constrain the Hilbcrt space. 
This can be done by examining the eigenvalue of a local 
operator rij. 

The J term is a reflection of the frustrated kinetic en- 
ergy at the level of the Hubbard model |}6) in the U — » oo 
limit, but at the level of the t-J model, the J term cannot 
be properly defined to be kinetic energy: it does not rep- 
resent motion of the particles described by the fermion 
operators. Moreover, it is neutral under gauge transfor- 
mation, because both Si and n, are. In contrast, the 
t-term is the kinetic energy; it picks up a Peierls phase 
under a gauge transformation and the constraint, being 
local, remains unchanged. Thus, it is meaningful to ask 
which term plays a more important role if the supercon- 
ductivity is described by the t-J model, but it is pure 
semantics to try to pin the mechanism down as being 
driven by kinetic as opposed to potential energy. What 
is potential energy at one level is kinetic at the other. If 
the t-J model is not adequate to describe superconduc- 
tivity, we must return to the Hubbard model, and the 
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partitioning of the kinetic and potential energies will be 
different. 

It is useful to examine the BCS theory of superconduc- 
tivity for which the effective hamiltonian is the reduced 
hamiltonian. A textbook calculation shows that the ki- 
netic energy is increased in the superconducting state, 
8(KE) = (A 2 /V)[l - N(0)V/2], while the potential en- 
ergy is lowered, S(PE) = —A 2 /V, due to the attraction 
of electrons mediated by phonons. Here A is the super- 
conducting gap, V the magnitude of the attractive in- 
teraction, and N(0) is the density of states at the Fermi 
energy. Although the phonon exchange is a kinetic pro- 
cess, its effect is correctly described as a potential energy 
at the level of the reduced hamiltonian. The increase 
of the kinetic energy is not in the least surprising be- 
cause BCS superconductivity develops on top of a Fermi 
liquid in which the kinetic energy is diagonal and un- 
frustrated. Therefore, it must necessarily be increased 
in the superconducting state. An interesting corollary 
is that if superconductivity is due to the lowering of the 
electronic kinetic energy in a suitable low energy effective 
hamiltonian, it could not develop on top of a Fermi liquid 
state, for in a Fermi liquid the kinetic energy operator is 
diagonal; the normal state will have to be a non-Fermi 
liquid. We may have a new class of superconductors, but 
it is still semantics to say that it is driven by kinetic en- 
ergy, for it will surely depend on the low energy effective 
hamiltonian, in which a part can appear as a potential 
energy, which could be a reflection of frustrated kinetic 
energy at the preceding level. 
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